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Abstract 

Oh! 

^ ' We present a new analysis of parton distributions of the proton. This incorporates a 

wide range of new data, an improved treatment of heavy flavours and a re-examination 
of prompt photon production. The new set (MRST) shows systematic differences from 
previous sets of partons which can be identified with particular features of the new data 
and with improvements in the analysis. We also investigate the sensitivities of the results 
to (i) the uncertainty in the determination of the gluon at large x, (ii) the value of 05 (M^) 
and (iii) the minimum cut on the data that are included in the global fit. 
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1. Introduction 

The last few years have seen a spectacular improvement in the precision and in the kinematic 
range of the experimental measurements of deep inelastic and related hard scattering processes. 
As a consequence the parton distributions of the proton are much better known, with tight 
constraints on the gluon and the quark sea for Bjorken x as low as 10~^. However, several 
significant sets of new data have become available recently which, when incorporated in a 
global analysis, will increase considerably our knowledge of the parton distributions. 

First let us summarize the new experimental measurements that have become available, and 
their implications. The new information includes the following. 

(i) We have new, more precise, measurements of the structure function F2' for deep inelastic 
electron-proton scattering by the HI and ZEUS collaborations at HERA ^]. The data 
now extend over a much wider kinematic range. Loosely speaking F2 and dF2/dlnQ^ 
serve to constrain the sea quark and gluonQ distributions in the region 10~^ ^ a; < 10^^. 

(ii) A re-analysis of the CCFR neutrino data has led to a new set of F^^ and xF^^ struc- 
ture function measurements 0. Besides affecting the quark densities, the new structure 
functions give a value of the strong coupling as much more in line with the world average 
determination than the original CCFR data 0]. 

(iii) There now exist measurements of the charm component of F2^ in electron-proton deep 
inelastic scattering at HERA 0, ^. These new data sample x ~ 10~^ and complement 
the existing EMC charm data with x ~ 10~^ Such data constrain both the charm sea 
and gluon distributions via the subprocesses 7*c — c and 'y*g cc respectively. 

(iv) Very precise measurements of prompt photon production, pp 7X, have become avail- 
able from the E706 collaboration . These data motivate us to look again at our treat- 
ment of this reaction, and in particular of the WA70 prompt photon measurements 0. 
We emphasize that such prompt photon data are the main constraints on the gluon[| out- 
side the HERA small x (x ~ 10^'^) domain, apart of course from the global momentum 
sum rule constraint. 



The E866 collaboration |T0[ have measured the asymmetry in Drell-Yan production in 
pp and pn colhsions over an extended x range, 0.03 ^ x < 0.35. The asymmetry data 
provide direct information on the x dependence of the difference, u — d, of the sea quark 



densities. Previously there existed only the single measurement of NA51 [jTT| at x = 0.18 



which revealed that d ~ 2m at this x value. Now much more information on m — d is 
available^. 



^At the lower values of there is a significant correlation between g{x, Q^) and the value of the strong 
coupling as(Q^). 

^Other reactions can in the future offer important constraints on the gluon, see (ix) below. 

^See also Section 7.2 where information on u — d from semi-inclusive deep inelastic data is discussed. 
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(vi) The CDF collaboration [12] have been able to considerably improve the precision of their 
measurements of the asymmetry of the rapidity distributions of the charged lepton from 
W"^ l^v decays at the Tevatron pp collider. The new data extend to larger values of 
lepton rapidity. These data offer a tight constraint on the u/d ratio of parton densities. 

(vii) NMC have now completed their structure function analysis, and supplied further data on 
^2*^, Fg^*^ and the ratio F!^'^ /F2^ The dedicated measurement of the ratio provides a 
valuable constraint on the u/d ratio. 



(viii) The data on Drell-Yan production obtained by the E772 collaboration cover a wider 



range of xp than the E605 data [jTSl which we have used to constrain the sea. For x^? ~ 



both experiments provide a useful measure of the quark sea at larger values of x, typically 
X ^ 0.3. For larger xp the E772 data probe, in principle, the valence quarks at x ~ 0.5 
and the sea quarks at x ~ 0.025. 

(ix) There are several other data sets which have the potential to provide important informa- 
tion on partons in the future. These include jet, W+ jet and heavy quark (6, t) production 
at Fermilab, and diffractive J/ip and dijet production at HERA, as well as jet production 
in deep inelastic scattering. 



The paper is organized as follows. In Section 2 we describe our procedure and the input 
parametrization of the partons. The optimum global set of partons, which we denote simply 
MRST, is presented together with two other sets which represent the possible range of be- 
haviour of the gluonQ. The quality of the description of the deep inelastic scattering data is 
shown in Section 3. The sensitivity to the cuts imposed on the data that are fitted and to the 
value of as are also discussed. In Section 4 we study the impact of the prompt photon data 
on the determination of the gluon. We pay particular attention to the transverse momentum 
(/ct) smearing of the incoming partons. In Section 5 we present a self-contained discussion of 
our new treatment of the heavy flavour (c, h) contributions to the structure functions. The de- 
scription of the data for the Drell-Yan process is given in Section 6 and in Section 7 we discuss 
the (m, d, s) flavour decomposition of the sea. Section 8 is devoted to the constraint imposed by 
the new W asymmetry data. Some implications of the impact of our new partons for processes 
observed at the Fermilab Tevatron collider are presented in Section 9. Finally in Section 10 we 
summarize the key features of our analysis. 



2. Global parton analysis 

The parton distributions fi are determined from a global fit to a wide range of deep inelastic 
and related hard scattering data. The basic procedure is to parametrize the /j(x,Q^) at a low 
value of = Ql such that the /j(x, Q"^) can be calculated at higher by using next-to-leading- 
order (NLO) DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) evolution equations. Data 

■^The FORTRAN code for all the parton sets described in this paper can be obtained from 
^ittp : //durpdg . dur . ac . uk/HEPDATA/PDF , or by contacting W.J.Stirling@durhain.ac.uk. 
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are fitted for all > Q\, where Q\ > Qq is a value of where perturbative QCD is believed 
to be the dominant contribution. We shall study the sensitivity of the results to variation of 
the choice of Q\. 

We parametrize the starting parton distributions at = Qq = 1 GeV^ where the number 
of active flavours is = 3. We work in the MS renormalization scheme and use the starting 
parametric forms 

xu, = A^x^^l - 3;)"^! + e„v^ + 7„a;) (1) 
xd, = AdX^^-'il-x^'il + edV^ + ^/dx) (2) 
xS = Asx-^^{l~xy^{l + esv^ + isx) (3) 
xg = AgX-^\l~x)'''{l + egV^ + -fgx). (4) 

The flavour structure of the light quark sea at Ql is taken to be 

2u = 0.45 -A (5) 
2d = 0.45 + A (6) 
2s = 0.2S (7) 

where the 2s/{u + d) ratio of 0.5 is chosen to obtain a s density consistent with the CCFR data 
on dimuon production see Fig. ^of Section 7. The parametric form of A, which specifies 



the difference between u and d, is taken to be 

xA = x{d-u) = A^x'^^il - xy'-^^{l + -fAX + 5ax'^). (8) 

The data require the integral Jq dx{d—u), which occurs in the Gottfried sum rule, to be positive. 
For the parton sets obtained in this analysis the integral is approximately equal to 0.1 giving 
a Gottfried sum 

/gs ^ - [Fi - ^ 0.27 (9) 



/o X 

12 ^ c n^\r2 



in the region Q ~ 5 GeV . The relevant data are the measurements of the asymmetry in Drell- 
Yan production in pp and pn collisions. The pioneering experiment of the NA51 collaboration 
111 measured the asymmetry at one value of x, x = 0.18. Very recently the E866 collaboration 



rUI have made measurements over a range of x. The new data indicate that, while indeed 
xA > for a; ^ 0.2, for larger x values A = d — u becomes small and may even become 
negative. The structure of xA is accommodated in our parametrizations by negative values of 
Sa- To ensure that the individual densities u and d are both positive at all x, we suppress the 
difference A at very large x by an extra factor of 2 in the exponent of (1 — x), see (8). 

For the first time our treatment of the heavy flavour densities, charm and bottom, is on 
a firm theoretical footing. These densities are determined by the other parton distributions 
and no extra parameters are introduced apart from the heavy quark masses. At very low 
the structure functions F^{x, Q"^), with H = c,b, are described by boson-gluon fusion and the 
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heavy quark densities turn on at ~ m|^. The procedure for ensuring a smooth continuation 
in the behaviour of in the threshold region is described in Section 5. 

A new feature of our analysis is the particular attention to the uncertainties in the gluon 
distribution at large x. The main constraints in this region are data on prompt photon produc- 
tion in or "pA collisions from the WA70 and the E706 experiments. The latter data 
confirm the implication from other high energy prompt photon experiments fl^ that a signif- 



icant initial state partonic transverse momenta is needed to obtain agreement with the NLO 
QCD prediction [0. This naturally raises the question whether such a transverse momentum 
component should be included when determining the behaviour of the gluon at large x{x ~ 0.4) 
from the lower energy (a/s = 23 GeV) WA70 measurements — even though the WA70 data 
can be adequately described without such a component]^. We findQ that the E706 data, which 
correspond to = 31.5 and 38.8 GeV, require the average value of the transverse momentum 
of the initial partonic system (/ct) ~ 1 GeV, and we expect this to be less for experiments 
at lower energies. We therefore begin by taking a canonical value of (A;t) = 0.4 GeV for the 
analysis of the WA70 data at -^i = 23 GeV. We then explore a range of gluon distributions 
which result from global analyses in which (/cy) goes from one extreme of (/ct) = to the 
other (/ct) = 0.64 GeV, which is the maximum value that we find compatible with a reasonable 
description of the WA70 data, see Section 4. We call the gluon distributions which correspond 
to (A;^) = 0, 0.4 and 0.64 GeV the higher, central and lower (large x) gluons respectively — 
since a smaller gluon density is compensated by a larger (/cr)- We denote the corresponding 
parton sets by MRST(5f |), MRST and MRST((7 |). Of course the inclusion of intrinsic (fcy) in 
this way can only be a crude approximation to the actual underlying physics, which presumably 
incorporates multiple emission of soft gluons as well as genuine non-perturbative higher-twist 
effects. However, unlike for the Drell-Yan process, there is as yet no complete theoretical treat- 
ment of this process, although a quantitative estimate of the effects of soft gluon emission has 
recently been made in Ref. ||I9|. However, the three choices of the large x gluon behaviour do 
give, we believe, a realistic indication of the uncertainty of the gluon distribution due to these 
effects. Table I lists the values of the parameters of (l)-(4),(8) for the three parton sets. 

The optimum global MRST description has the QCD parameter Aj^(n/ = 4) = 300 MeV, 
which corresponds to asiM'z) = 0.1175, in excellent agreement with the world average value 
as{Mz) = 0.118 [0. The same as value is also used for the MRST(5f | and g |) partons[]. 
In the literature errors between ±0.005 and ±0.003 are quoted on the world average value of 
as{M^). In Section 3.3 we will present parton sets which cover the range given by the more 
conservative error of ±0.005, in each case taking (/c^) = 0.4 GeV when analysing the WA70 
data. 

From Table I we see that the values of are negative, which imply a 'valence' type of 
behaviour for the gluon at low x for = 1 GeV^. Care should be taken not to attach physical 

^In previous MRS analyses of the WA70 data we did not include initial state transverse momentum. 
^The calculation is described in Section 4 and compared with existing descriptions of the E706 data. 
''The optimum value of as{M^) for the (g t and g [) parton sets is close to that of the central gluon fit. 
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significance to this behaviour, as it arises from an 'extrapolation' outside the domain of the fitted 
data. Indeed as increases the behaviour rapidly changes so that by = 2 GeV^ the gluon 
distributions are approximately 'flat' in x. Fig. [I| shows the three gluon solutions as varies 
from 2 to 100 GeV^ and we can see how evolution up in soon blurs the distinction between 
the initial starting distributions. Notice also that at low the distinct behaviour at a; ~ 0.4 
is compensated by the opposite behaviour at a; ~ 0.05 so that in each case the momentum 
fraction carried by the gluon remains at roughly 35% at Qq = 1 GeV^. This structure is more 
evident in Fig. |^ which shows the ratio of the gluons at = 10 and 10^ GeV^. We see that 
all our three gluon distributions converge for x ^ 0.01 due to the requirement of fitting the 
HERA data. The CTEQ4M ||21| gluon distribution is also shownQ and the comparison will be 
discussed later. 

Table II lists the fraction of momentum carried by the individual partons as a function of 
for the central solution, MRST. For = 200 GeV^, for example, 46%, 31% and 23% is 

carried by the gluon, valence and sea quarks respectively. In fact the flavour decomposition of 

the sea momentum fraction is 

u, d, s, c, b = 7, 8, 5, 3, 1%, (10) 

demonstrating the growth of the strange, charm and bottom distributions with increasing Q^. 

The wide range of processes used in the global analyses is listed in Table III. We include deep 
inelastic scattering (DIS) data from HI 0], ZEUS§, BCDMS[g3, NMCfg, E665l2|, SLAC12§ 



and CCFR§. In general we only fit to DIS data with Q'^ > 2 GeV^ and VV^ > 10 GeV^, but 
in order to include very small x measurements of we admit the HERA data for down to 
1.5 GeV^. Compared to our 1996 global analysis |2^, the HERA data are updated to include 



the HI shifted vertex as well as nominal vertex data, the NMC data now include all five beam 
energies and, finally, we use the reanalysed CCFR neutrino data. Table IV shows the values 
for all these DIS data for the three 'gluon' fits described above. The important constraints 
from non-DIS data are discussed in detail in later sections. Figures ^ and ^ show the MRST 
parton distributions as a function of x for = 20 GeV^ and = 10^ GeV^ respectively, 
while Fig. |^ compares them with those of MRS(R2), our favoured set of the previous analysis 
We discuss the comparison of the new partons with previous sets in Section 10. 



3. Description of DIS data 

The description of the DIS structure function data by the MRST partons is shown in Figs. ^ 
1^. Overall the quality of the fit is satisfactory as reflected by the values listed in Table 
IV. However, it is informative to note special features of the fit and, in particular, to highlight 
those areas in which the description of the DIS data is systematically poorer than average. 

A comparison with the small x data that are used in the fit is shown in Fig. ^. For the 
purposes of illustration, data at adjacent x values are grouped together at a mean value, together 

^Very recently CTEQ have attempted to estimate the uncertainty on the gluon, without using any 
constraints from prompt photon data. 
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with the MRST fit. Very recently the HI collaboration [27| have made available unpublished 
preliminary measurements of F2 from the 1995 and 1996 runs. These data are not included in 
the fit but a comparison with MRST is shown in Fig. ^ It is clear that the data in Fig. ^ are 
sufficient to put strong constraints on the small x behaviour of both the sea quark and gluon 
distributions. Loosely speaking F2 and dF2/ dhiQ"^ determine the exponents, A5 and Ag, of 
the sea and gluon distributions respectively, as well as constraining the overall normalization. 
The values of A5 and, particularly, A^ depend sensitively on the value of = Ql chosen to 
parametrize the input distributions. 

We may look at how the exponents A5 and Xg vary with by fitting the sea quark and 
gluon distributions of the MRST partons to the forms 

xMx,Q') = A(Q2)x-^'(«') (11) 



as X ^ 0. The results are shown in Fig. 13. As increases from the input scale, Qg = 1 GeV^, 



we see that the valence-like character of the gluon rapidly disappears due to evolution being 
driven by the much steeper sea, and that by — 2 GeV^ the gluon is 'flat' in x, that is A^ = 0. 
By — 6 GeV^ we see that A^ = A^, which incidentally is close to the assumption]^ made in 
the early NLO global analyses in which the input scale was chosen to be Qg = 4 GeV^. For 
higher values of the gluon exponent 'leads' that of the sea, A^ > A5, since the gluon drives 
the sea quark distribution via the g —>■ qq transition. 

Fig. is an alternative way of looking at the quality of the description of dF2/ d\n.Q^ at 
low X. The continuous curves on the plots show the values of the slope dF2/d In versus x for 
both the HI and ZEUS data, compared with the slope found in the MRST fit (evaluated at the 
particular values of appropriate to the experimental data). Though the overall description 
is satisfactory, it is possible that for x ^ 10~^ there may be a systematic difference between the 
data and the fit which reflects the onset of In 1/x contributions which are outside the scope of 
our NLO DGLAP analysis. This systematic trend is even more evident in the preliminary HI 
data from the 1995/96 runs 

13- As the precision of the HERA measurements of F2 improves, it 



will be interesting to see whether or not the statistical significance of the discrepancy increases. 

The description of the NMC DIS data for F| and F!^ is shown in Fig. | and Fig. |. It 
is apparent that the data have systematically a larger slope, dF2/d\nQ'^, than the fit. This 
is well illustrated by the continuous curve in Fig. ^ which shows dF2/dlnQ^ for the NMC 
F2 data. The discrepancy indicates that the NMC data would prefer a larger gluon in this x 
region and/or a larger as value than that of MRST. In turn the larger gluon would imply (by 
momentum conservation) a smaller gluon in the very small x domain, contrary to the HERA 
data. The fit is a compromise between these data sets, but it demonstrates the tight constraints 
now imposed by the increased precision of the F2 data. A particular virtue of the NMC data 
is the accurate measurement of F2/F2. We postpone a study of their implications until we 

^Approaches in which Ag is tied to at lower scales (say Qg < 2 GeV^), as in the dynamical GRV model 
, will clearly have difficulty in fitting the new HERA data. 
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discuss the description of the asymmetry in Drell-Yan production in pp and pn colhsions and 
the rapidity asymmetry in the processes pp — W^X. 



The BCDMS data [23| cover the range 0.1 < a; < 0.75 and are the most precise data at large 
X, see Fig. |1^. From the figure it is apparent that the data would prefer a smaller value of 
as than that found in the global fit, namely ^^(Ml) = 0.1175. Indeed if the BCDMS data 



are analysed on their own (apart from including SLAC data |25| to constrain the higher twist 
contribution) then a5(M|) is found to be 0.113 ± 0.005 [^]. Our optimum value of as is 
therefore within a standard deviation of the BCDMS determination. 

The re-analysed CCFR neutrino measurements of ^2^^ and xF^^ are compared with the 
MRST values in Figs. |Tl] and [1^. The long-standing discrepancy between the CCFR and 
the NMC Fg''' measurements for x < 0.1 remains: in this x region the measurements are 
in excess of the Fg^'^ data by a significantly larger amount than that implied by the strange 
quark distribution]^ determined from dimuon data, see Fig. ^ of Section 7.1. However, it is 



important to note that the neutrino data must be corrected for heavy target effects. In Figs. |TT 



and |T2| we have subjected the MRST curves to a heavy target correction. The parametric form 
that we use for the heavy target correction factor i?HT is deduced from a Q^-independent fit to 
the EMC effect for the scattering of muons on a heavy nuclear target {A = 56). Explicitly, 

1.238 + 0.203 logio ^ for x < 0.0903 
Rut = { 1.026 for 0.0903 < a; < 0.234 (12) 

0.783- 0.385 logioa^ for 0.234 < x. 

Note that the correction factor that we obtain in this way is more severe at low x than that 
implied by shadowing. This is one reason why the (dashed) curves are considerably below the 
X < 0.1 neutrino data. It is not clear whether the correction factor should be the same for 
neutral current and charged current DIS data, or be the same for F2 and xF^ neutrino data. 
For these reasons we do not include the CCFR heavy target data for x < 0.1 in the fit. The 
remaining neutrino data are well described. However, it is possible to see that the fit slightly 
underestimates the slopes as a function of InQ^, which refiects the value a5(M|) = 0.119 
obtained from fitting to the neutrino data alone 0. 

3.1 Fl and implications for partons 

In the DIS experiments it is not the structure function F2{x, Q^) which is measured directly 
but the differential cross section. Defining the rescaled differential cross section 

" 2^a2 [l + (l-y)2]dxdQ^' ^^^^ 

where y = Q^/xs, we have 

a(x, Q2) = F2ix, Q^) - f Fl{x, Q'). (14) 

[1 + {l-yy\ 



lORecall, at LO, that xs{x) = ^F^^{x) - 3F^'^{x). 
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Since both y and are usually small the latter term in this expression is usually negligible, 
and the measurement of F2 is effectively direct. 

However, the analysis of data on the longitudinal structure function Q^) is in principle 
an important probe of the parton distributions. This is particularly the case for the gluon at 
small X since in this region, to a good approximation, we have the relationship 

xg{x, Q^) = 5.9(3a5/47r)Fi(0.4a;, Q^) (15) 



for three massless flavours at leading order in as (though the next-to-leading correction 



leads to an increased gluon relative to a fixed Fl). Nevertheless, until recently there have 
been little data on Fl, and these have been at high x (where there are likely to be important 
higher twist effects) and have very large errors. The situation is now beginning to change. The 
HERA experiments measure the differential cross section a.t y > 0.5, and are thus sensitive 
to the component due to Fl. Also, the NMCollaboration have direct measurements of Fl for 
0.1 > x > 0.01, obtained by data runs for different beam energies, and hence different values 
of y. 

The only consistent way in which to analyse the HERA data at large y is to calculate both 
the NLO expressions for F2 and Fl (i.e. using the C(a|) coefficient functions for Fl P^, 153 



where those for heavy quarks use the prescription of ||4g] ) and to compare ([1^) with the measured 
a. Of course, this is equivalent to the correction of the extracted values of F2{x,Q'^) to take 
account of the predicted values of Fl{x,Q'^), and it is this latter procedure that we employ 
when fitting to the HERA data. (This results in corrections of at most 2 — 3% to the values of 
F2{x,Q'^) quoted in [||, Q.) Since in the pubhshed data the value of y is nearly always < 0.6, 
there is relatively little sensitivity to the value of Fl- However, the HI collaboration have also 
published a number of measurements of a for y = 0.7 [Q, and in the preliminary 1995/96 data 



have reproduced these measurements, and also produced measurements at y = 0.82 In 
Fig. ^ we show a comparison of the prediction for a{x, Q"^) obtained from the MRST partons 
with these preliminary HI data. It is clear that the y = 0.82 points (the first data point in 
each plot) lie below the theoretical curves in general (which implies that the predicted Fl is 
too small). However, it is also clear that the points for 0.01 ^ x ^ 0.001 tend to lie above the 



curve. Indeed, the curves in Fig. |16| resemble the curves in Fig. |1J: there is a tendency for 
dF2{x,Q'^)/dhiQ'^ to be too small at x ~ 0.005 and hence F2{x,Q'^) tends to be too small at 
high Q^, and at a; ~ 0.0005 there is a tendency for dF2{x, Q^)/d\iaQ'^ to be too large and hence 
F2{x, Q^) tends to be too large at high Q^. Thus, also bearing in mind the large errors on the 
points at y = 0.82, we feel that it is premature to claim any inconsistency in the NLO prediction 
for Fl{x, Q^)-f^ Alternative theoretical treatments to ours, in particular the inclusion of leading 
ln(l/x) terms, lead to different predictions for Fl for similar fits to F2, and measurements at 

^^We note that in the NLO fit performed by HI there is no direct constraint on the gluon at large x. From 
the momentum sum rule this leads to more flexibility for the gluon at small x, and hence the details of their 
best fit are rather different from ours in this region. 
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high y are therefore an important test of such approaches.0 However, direct measurements of 
Fl would provide an even better test. We exhibit the predictions for in the HERA kinematic 
range obtained using the MRST partons in Fig. 

In Fig. |l^ we compare our predictions for Fl{x,Q'^), made using the partons resulting 



from each of the three parton sets, MRST, MRST(5f |) and MRST((7 J,), with the direct mea- 
surements made by NMC||13[. There is little variation in the predictions, and each provides a 
perfectly satisfactory description of the data. Hence, the NLO calculation of structure functions 
seems to be compatible with both direct and indirect data on F^. 

3.2 Sensitivity to cuts on the data fitted 

Given that there are potentially important higher twist and ln(l/x) contributions at small 
and/or small x it is instructive to explore the sensitivity of our fits to the minimum 
and/or x cuts. Recall that our minimum cut is a.t Qf = 2 GeV^, except for the HERA 
data where the cut is at 1.5 GeV^. We have imposed no minimum cut on x. We have made 
repeated global analyses with different minimum cuts up to a value Qf = 10 GeV^. Note 
that increasing Ql has the effect of removing much of the lowest x data from the fit. The dashed 
curve in Fig. shows the effect on dF2/dlnQ^ when the cut is increased to = 10 GeV^. Not 
surprisingly the description of the slope determined from the HERA data with > 10 GeV^ 
is much improved. We should note that with Qf = 10 GeV^ a significant fraction of the NMC 
data is also excluded, see Fig. This removes a strong constraint on the behaviour of the 
gluon at intermediate x thus allowing the gluon to increase at small x, which improves the fit 
to the HERA data. As we can see from Fig. |I9] the gluon obtained from the Qf = 10 GeV^ fit 



is larger at small x than the standard MRST gluon. As expected the difference between the 
gluons decreases rapidly with increasing Q^. The effect of an intermediate choice of Qf can be 
anticipated by interpolating the Qf = 2 and Qf = 10 GeV^ results. 

Fig. |2y shows the sensitivity of the parton distributions at = 10 GeV^ to the minimum 
cut on the structure function data that are included in the fit. Recall that only data with 



Q > Ql are included. Fig 20 compares the results of analyses with cuts at Q| = 5 and 



Ql = 10 GeV^ to our defauh MRST set obtained by taking the cut at = 2 GeV^. The plot 
is interesting because the variation of the values of the partons with the Qf cut reflects the 
interplay of the constraints imposed by the various data sets. For instance the u and d quarks 
for 0.01 ^ X < 0.5 are stable to a choice of Ql in the range 2-10 GeV^, since precise data remain 
in this domain even for the highest Ql cut. For the higher Ql values the small enhancement 
(up to at most 2%) in the u distribution in the region 0.01 < x < 0.1 can be understood by 
looking at the description of the highest Q^ NMC F2 data points in this x range, see Fig. ||. 
The variation of the gluon with Ql has been discussed above, and is responsible for the similar 
variation of the charm distribution. The decrease of the light quarks at small x [x < few x 10~^) 
with increasing Ql is partly due to compensation for the increase of the charm contribution to 



^^For a discussion of such tests, and in particular a demonstration that the high y data gives strong evidence 
against the vahdity of a LO-in-ag fit, see ]35[. 
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and partly due to the decrease of itself induced by the larger slope dF2/ dln.Q'^ required 
by the higher data. 

We have also repeated the global analysis with various x cuts on the data up to x^^^i = 0.01. 
This removes more than half of the HERA data, but leaves the NMC and other fixed target 
data virtually untouched. The removal of the constraint on the gluon at very small x allows a 
larger gluon in the region x ~ 0.1, leading to an improvement in the description of the NMC 



data illustrated in Fig. Nevertheless, as can be seen, the improvement is not as great as 
might be expected. 

3.3 Sensitivity to as 

We have mentioned that as{M^) = 0.1175 yields the optimum of the global fit to the 
combined data sets. Next we explore the sensitivity to the variation of the value of a^. To do 
this we perform global analyses with fixed values of as in the range ±0.005 of our optimum 
value. In each case we use (fcr) = 0.4 GeV when analysing the WA70 data, which corresponds 
to the central gluon distribution of the three 'gluon' fits described in the previous section. The 
contributions to the total coming from the various data sets are plotted as a function of as 



in Fig. 21 . We emphasize that this is not the optimum for a particular data set fitted on its 
own, but is the contribution to for the global fit which, of necessity, has to make compromises 
between the descriptions of the various data sets. As expected from our previous discussion, we 
see the opposite trend for the of the BCDMS data (which favour a smaller as) and the CCFR 
data (which favour a larger as)- Similarly the NMC data favour a larger as to compensate for 
the global fit yielding a smaller gluon than that which would be obtained by fitting to the data 
set on its own. It is noticeable that the recent, more precise, HERA data give a contribution to 
which is less sensitive to variation in as than the earlier HERA measurements (see ref. [p6| ). 
In summary we see that the overall minimum value as{M'^) = 0.1175 in the global fit is a pinch 
between the BCDMS and HERA data favouring smaller as values and the NMC, SLAC and 
CCFR data favouring larger values. 

An independent sensitive measure of the value of is the single jet inclusive Ex distribu- 



tion measured in the Fermilab pp experiments ||3^, We shall see in Section 9 that these 
data favour values of as'(Ml) in the region 0.115-0.118. There thus seems to be general agree- 
ment that the value of a5(M|) is in the region of our optimal value 0.1175 with a spread of 
about ±0.0025. It is useful to have a range of parton sets available for different values of as, 
so we present four additional sets which cover a conservative range of ±0.005 about our opti- 
mal value. We denote these by MRST(a5 ij), MRST(a5 i), MRST(a5 T) and MRST(as IT) 
corresponding to as{M^) = 0.1125, 0.1150, 0.1200 and 0.1225 respectively. 

Given the recent interest in the very high region at HERA it is important to quantify 
the uncertainty in the extrapolation of Ff' to high at a large value of x. In this x region the 
main effect comes from the uncertainty in the value of ^5 in the DGLAP evolution of Ff'. The 



effect is illustrated at x = 0.45 in Fig. which shows the spread in the extrapolated values 
of F2 arising from as varying across the interval 0.1175 ± 0.005. 
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4. Prompt photon production and the gluon at large x 

In previous MRS parton analyses the WA70 data for pp — ^ 7X have been a key constraint 
on the gluon distribution for x ~ 0.3-0.5. In these analyses we have not included any initial 
state partonic transverse momenta, that is we have taken (fc-r) = 0, in fitting to the prompt 
photon data. If we continue to fit the WA70 photon pt spectrum in this way then the global 
analysis yields the set of partons that we have called MRST(5f |) in Table I. We use the MS 
renormalization and factorization prescriptions, with a common scale Q = Pt/2. We perform 
a full NLO calculation including the effects of fragmentation A good description of the 

WA70 data is achieved. 

However, there is now compelling evidence for the need to include non-zero transverse mo- 
mentum kr of the incoming partons (arising from parton multigluon emission^ and from non- 
perturbative 'intrinsic' partonic transverse momentum). The reason is apparent if we consider 
all the data for prompt photon production in high energy pp and pp collisions simultaneously. 
Collectively these data span the entire interval 0.1 ^ x ^ 0.5. A major experimental challenge 
in these experiments is to cleanly extract the prompt photon signal from the copious back- 
ground of 7r° and rj decays. There is a pattern of deviation between theory and experiment in 
the shape of the photon px spectrum. The data fall off more steeply with increasing pt than 
the NLO QCD predictions. Neither changes of scale nor the introduction of fragmentation 
effects can resolve the discrepancyf^ since the various experiments probe different ranges of 



X xt = 2pT I y/s. On the other hand it has been shown ^l| that the discrepancy can be 
removed by a broadening of the initial state parton kx which increases with the energy y/s. As 
already mentioned in Section 2, there is no complete theoretical treatment currently available 
which would allow a parameter free description of the broadening of the pt distribution. There 
is evidence from vr^vr'^ and 77 production that an approximately Gaussian smearing form repro- 
duces the observed broadening (see, for example, Ref. Q]) and that the width of the Gaussian 
increases with energy. 

The parton kx effect is found to be least in fitting the data due to WA70 — the lowest energy 
prompt photon experiment. Moreover these data do not exhibit the px shape discrepancy with 
QCD that is seen in the other experiments [0, § . For these reasons we have set (/cr) = in our 
previous analyses. However, if we include this type of analysis in our new global fit and use the 
resulting partons, MRST(5f |), to predict the high precision E706 prompt photon pt spectra 
the description is disastrous. To reconcile our prediction with the E706 data we may fold in a 
Gaussian kx spectrum^ with {kr) — 1 GeV, or to be precise {kx) = 0.87(0.97) GeV for data 

^^We thank Werner Vogelsang for performing the relevant calculations. 

^^An estimate of the amount of smearing expected from the perturbative component based on the resummation 
of leading logs has been made in Ref. jlOj . 

i^Vogelsang and Vogt |Q have demonstrated that these effects can improve the description of a single 
prompt photon experiment. However, experiments at different ^/s reproduce a similar pattern, but in different 
X intervals. 

-'^^This is in qualitative agreement with the findings of the E706 collaboration |Q, although the detailed 
prescriptions for the kr smearing are different. We smear the perturbative QCD distribution by first making an 
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taken at laboratory momentum piab = 530(800) GeV. As a consequence it can be argued that 
our (fc-r) = analysis of the WA70 data is inconsistent. We should include {kx) 7^ in the 
description of these lower energy data but with a smaller value|^ of {kj) than that needed to 
describe the E706 data. 

We therefore repeat the global analysis but fit to the WA70 data using a Gaussian partonic 
/ct spectrum with (/ct) = 0.4 GeV, corresponding to 280 MeV per incoming parton. The 



description of the WA70 data is shown by the continuous curve in Fig. ^ For comparison the 
dashed curve shows the unsmeared prediction, which of necessity undershoots the data. The 
resulting set of partons are labelled simply MRST in Table I, and give an equally good fit to 
those of MRST((7 f) with {kx) = 0. How well are the E706 photon data described by the MRST 
partons, and in particular by the new gluon which is smaller at x ~ 0.4? From the continuous 



curves in Figs. |2^ and |^ we see an excellent description of the piab = 530 (800) GeV E706 data 
is obtained if we take {kr) = 0.92 (1.01) GeV. 

We can regard the MRST(5f |) partons with (kx) = as one extremum. Conversely how 
large can we take (fcy) to be to describe the WA70 data and still retain a satisfactory fit? We 
find that we cannot choose (kx) to be arbitrarily large because not only does the gluon become 
smaller, but it becomes steeper and eventually the shape of the pt spectrum is not reproduced. 
The maximum value of (/ct) for which a reasonable fit to the WA70 data can still be found is 
0.64 GeV. The parton set corresponding to this upper extremum for (fcr) is labelled MRST(5f |), 
to indicate that it has the smallest gluon at x ~ 0.4. The E706 data are well described by the 
MRST(^ i) partons provided we take (kr) = 0.97 (1.04) GeV at piab = 530 (800) GeV. 

So far we have considered the variation of the partons, and in particular of the gluon, due 
to the uncertainties in the (fcy) smearing. Our preferred set of partons with (/cy) = 0.4 GeV is 
MRST. The 'extremum' parton sets with (kx) = and 0.64 GeV are labelled MRST((yf |) and 
MRST(5f I) in Table I, with gluons which are respectively larger and smaller than the MRST 
gluon at X ~ 0.4. 

There is also a non-negligible dependence on the choice of scale. For instance the effect of 
changing the scale from Q = /2 to Q = pt is shown in Fig. We see that the unsmeared 
cross section is decreased by some 30%, which can be compensated by a relatively modest in- 
crease in the size of the gluon distribution and/or in (/ct)- That is the effect of the change of 
scale is considerably less than the uncertainties associated with smearing. In principle the for- 
mer can be reduced by a knowledge of the NNLO perturbative contributions, while a reduction 
in the latter will require a more detailed theoretical understanding of the origin of the partonic 
transverse momentum. 



analytic continuation for px < pa = 3 GeV of (dcr/dp^)QCD of the form exp(^ -^q aip^j.) to regulate the infra-red 
singularity at px = 0, and then we convolute with a Gaussian form (I/ttct) exp(— fcy/cr) where a = {4 / tt) (kx)'^ ■ 
The results are not sensitive to the particular choice of po- 

^^The empirical evidence that {kx) increases with ^/s is supported by a similar effect in Drell-Yan production 
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5. Treatment of Heavy Flavours 

Until recently the treatment of heavy quark distributions in MRS and CTEQ global analyses 
has been rather naive. In previous analyses the charm and bottom quarks were regarded as 
infinitely massive below a threshold = m^, and then being treated as massless, and thus 
evolving according to the normal massless evolution equations above this threshold. Up to NLO 
in as this prescription guarantees that the correct results will be obtained asymptotically, but 
is clearly rather unsatisfactory near threshold where there should be a smooth threshold at 
= Q'^{l/x — 1) = 4m|^, where W"^ is the invariant mass of the hadronic system, rather 
than an abrupt threshold in = rnjj. Nevertheless, choosing the slightly high value of 
= 2.7 GeV^, a reasonable match to the EMC data on the charm structure function for 
> 4 GeV^ was obtained^, and the low contribution due to charm for the total structure 
function rendered a more complete treatment of heavy quark contributions unnecessary in this 
X and range. 

An alternative procedure to that outlined above is where all charm is regarded as being 
produced from the hard scatter between the electroweak boson and a light parton, i.e. the 
number of active flavours, nj, is 3 and the charm cross section is generated (mainly) by photon- 
gluon fusion (PGF). This corresponds to the so-called fixed fiavour number scheme (FFNS) and 
it incorporates the correct threshold behaviour automatically. For example, at order 0{as) the 
charm structure function is given by 



F^{x,Q\ml) = ^^c«^^(gV^D ® gn^f^'), (16) 

ZTT 



as(jJ^ 
2n 

where the coefficient function (CF), which is convoluted with LO evolved gluon density gnf=3, 
is 



(Pi(--).4^-ii-3--)-8(;|)-) i„(i±i; 

2 



+ {8z{l - z) - 1 - A--^z{l - z))v 



e{W' - Ami) (17) 



where = Q^{l/z — l) is the gluon quark centre of mass energy, v is the velocity of the charm 
quark or antiquark in the photon-gluon centre-of-mass frame, defined by f ^ = 1 — Aml/W'^, 
and Pqg{z) = z^ + {1 — z)^, the LO quark-gluon splitting function. These t>-dependent terms 
ensure that the coefficient function tends to zero smoothly as W'^ = Ami is approached from 
below, and hence the structure function has a smooth threshold in W"^. This method does not 
sum potentially large logarithms in Q"^ /ml, and thus is unsuitable for Q"^ ^ m^, but provides 
an acceptable description provided is not large and one is not interested in the concept 
of a charm quark density. It is the method used to produce the most recent GRV structure 
functions |23], and is also used in the analyses by HI and ZEUS. 



^*In the MRS global analysis of |£6| the charm evolved from the low value of = 1 GeV^ but was suppressed 
by a phenomenological damping factor. 



13 



However, the more recent measurements of charm production at HERA ^, ^ emphasise 
the importance of having a consistent theoretical framework for heavy flavour production in 
deep inelastic scattering. Not only are there more direct measurements of the charm structure 
function ^2^^, but the charm contribution could be 20% or more of the total F2 at small x. 
Indeed, even the NMC data which contains only 5 — 10% charm, but has rather smaller errors 
than the HERA data, is sensitive to the treatment of charm. Hence, a modern global analysis 
of structure functions must necessarily include a satisfactory description of 



As a consequence there have been several recent theoretical studies p4| , |46[] to improve the 
treatment of heavy quark mass effects in deep inelastic scattering. For instance ref . |^ proposes 
a simple procedure to sum up the leading (and next-to-leading) log contributions of Feynman 
diagrams including explicitly the mn 7^ mass effects. It is straightforward to generalize this 
procedure to any order. In this approach the natural scale to resolve charm quarks in the 
proton is = (m^ -|- k^)/{z{l — z)) > 4m^, whereas the conventional MS scheme, which we 
adopt, requires the charm threshold to be at = m^. {kx and z define the momentum of the 
charm quark). For this reason the procedure is difficult to reconcile with the MS scheme. 

5.1. Theoretical procedure 

In order to have a reliable treatment of massive quarks over the whole range of we 
must clearly use an approach which extrapolates smoothly from the FFNS at low to the 
massless evolution at high Q^, maintaining the correct ordering in both schemes. To do this 
we use a method which has recently been developed by two of the authors and is discussed in 
detail in and more briefly in [0. Since this treatment of charm is such a major change 
to our previous analyses we present the method briefly here. First we note that in the FFNS 
(|T6|) is valid up to corrections of (9(A^/m^) while the massless prescription is valid only up 
to corrections of (9(m^//i^), i.e. threshold corrections. In order to improve the accuracy of 
the latter scheme we need to examine the connection between the parton densities in the two 
schemes. The connection between the MS parton densities for 3 and 4 flavours takes the form 

c+{z,^i^,fiyml) = A'^fiyml) ® (7„,=3(/i') 
gnf=4{z,i^\fiyml) = A^3{^i^/ml) ® ^„^=3(/i^) (18) 

at leading order, where the elements A^"- which contain ln(/i^/m^) terms, are, in general, part 
of a full 5x4 matrix which also connects the light quark flavours. Hence the charm distribution 
c+ = c -|- c is determined entirely in terms of the light parton distributions, and it is the above 
equations which lead to the requirement of evolving from zero charm at /i^ = m^. Since we use 
the scale choice /i^ = for all the light partons we also take this simple choice for the heavy 
quark structure function. Thus, from now on we will always use instead of /i^. 

For ^ the equivalence of the FFNS and the massless scheme at all orders lead to 
the connections between the CF's in the two schemes up to C(m^/Q^) [Q, in particular up to 
0{al) 

Cf{z,Q^lml) = C:f=' A^^iQ^ml) 

+ C;f=^ A^^Q^ml) + 0{mllQ^). (19) 
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The details of the connection are fully worked out in ||^. To improve the accuracy of ([T9|), where 



the uncertainty is reduced to C(A^/m^), requires defining 'corrected' CP's, C^^ (i = 1, . . . ,4), 
in another Uf = 4 scheme - the variable flavour number scheme (VFNS) - where one can write 



where the corrected CF's are related to the FFNS CF's by 



+ CYiQVml) ® guMQ'^Q'lm^) + 0{K'/ml), (20) 



Cr{z, Q'/ml) = CY{Q'/ml) ® A^\QVm% (21) 

the new nf=4 CF's now being exact at all values of Q^. 

Hence, our procedure is to use the FFNS for < ml where it should be very reliable and 
switch to the VFNS for > m^. (The precise choice of the transition point is undetermined, 
however, taking = removes complications arising from ln((5^/m^) terms in the matching 
conditions between the partons at threshold.) In order to define the VFNS one must solve (|2T|) 
for the . Unfortunately the all-orders matching of F2 in the two schemes, from which (|2T|) 
arose, is not sufficient since, for example, at low orders the single quantity Cj^ is expressed in 
terms of the two quantities C^^ and Cj^. We stress that any choice satisfying ( |21| ) is 'correct' 
in the sense that it leads to the same all orders expression. Nevertheless, each choice leads to 
a different expression if one uses the usual rules of combining coefficient functions and parton 
distributions of a given order to obtain a fixed order in as expression for the structure functions. 
In order to remove this ambiguity we apply a sensible, physically motivated constraint and 
impose not only continuity of the structure function but also demand, in addition, order-by- 
order matching of the evolution of F| at threshold. 

The exphcit form of (PI) at 0{as) is 



gnf=4{z,Q^,Q'^/ml) = gnj=3{z,Q'^) - ^ \n(^\ gnf=3- (22) 



Inserting the implied expressions for the matrix elements A'^^{z, Q"^ /ml) and A^^{z, Q"^ /ml) into 



([19|) gives the relation (first seen in 

C« ^^(., Q'/ml) = C« v^(., Q-/ml) + "''(Q'/ml) ® ln(g) (23) 

connecting the gluonic CF's in the FFNS and VFNS. Let us now consider the evolution of F2. 
From (|16D the LO expression in the FFNS for the In derivative is simply 

dh[Q^ 2^ rfW ^ ^ 
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The corresponding expression obtained by differentiating the LO expression in the VFNS, for 
just above m^, is 



rflng2 ^ 



Pi ® 9nr- 



pO 



c4Q2) . (25) 



At Q 



^ — m^, the terms in 



involving c+ vanish because of 
asition point imr 



and so demanding continuity 



of the evolution across the transition point immediately leads, from (^,25), to 



(26) 



(ilnQ2 

Generalising this relation to be the definition of Cf"* ^^(-z, Q"^ / m^) at all guarantees a smooth 
passage for charm structure function from < to > ml, by definition. It is also easy 
to see that in the limit Q"^ ^ oo, 



dC^'^''^{z,Qyml) 
rflnQ2 



(27) 



Hence, from (P^D, we see that C^°^ ^^(2;, Q^/m^) must indeed tend to the usual simple form 
z (5(1 — z) in this limit. Also, since C^^-* ^^{z, Q'^/ml) contains the factor 9{W'^ — 4m^) so does 
its In derivative, thus ensuring the correct threshold behaviour in for C^P^ and in turn 
for F2 at LO. Furthermore (^61) allows the gluonic CP in the VFNS to be written as 



(28) 



C^'^^^iz,Qyml) = Cl^^^^iz,Qyml) - '^^''"^[nQ^ ''"^'^^ ^"f^^ 



and g^jgQ j^g^g ^YiQ correct threshold behaviour as Q'^/ml 00 and C^^-' ^^(z, Q^/m^) 

tends to the correct asymptotic MS limit. The extension of this procedure to any arbitrary 
order, i.e. continuity of the derivative in the gluon sector, is described in full in [^ . 



The implementation of the charm coefficient function is also described in detail in [^, and 
results in the relatively straightforward expression 



^0 ^JCj^^)^^{z,Qyml) fxY dc4x/z,Q') 



+ 



dx 



dz 





dlnQ^ 


dCi^^ 






d\nQ^ 


dCi^) 




dlnQ^ 



d{x/z) 
c+{x/z,Q'^) 



dz'r{z') ^c+{x/zz',Q^ 

Jx/z ZZ 



where xq = (1 + Am^/Q ) and r{z) is given by 



r(z) 



1 

Z2 



K 3 . /V7 1 

cos — In - H -= sm — m - 

V 2 V7 V 2 z 



z' 
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(29) 



(30) 



This general method can be apphes at all orders, and the 0{as) charm coefficient function is 
determined by demanding continuity of dF^i^x, , ml) /dlnQ"^ (in the gluon sector) at (9(q;|), 



and is discussed in detail in However, in practice its contribution to the charm structure 
function is only at the level of a couple of percent at most, and it can be treated using a 
phenomenological approximate expression. 

Thus we can calculate the charm structure function at NLO. For < we use the usual 
FFNS expression, i.e. using coefficient functions |Q to C(a|)[^ and parton distributions with 
3 light quarks. For Q"^ > ml we use the VFNS coefficient functions to 0{as) and the partons 
are evolved via the NLO DGLAP equations in MS scheme with 4 massless quarks. Since the 
coefficient functions reduce to the usual massless expressions as Q^/m^ oo, the structure 
functions approach the previous massless expressions in this limit. 

We feel we should distinguish between this approach and previous implementations of a 
VFNS. In g (H) served as the definition for CW vf ^^^^^ pQp Qp ^-^^^ 

assumed form of Cf^ given by 

^^(., Q'/ml) = z Sixo - ^) (l + ' ^0 = (l + ^) (31) 

where the delta-function describes the tree-level diagram for a massive quark scattering from 
a photon and the modified argument of the delta-function follows from demanding that the 
massive quark is on-shell.[^ We believe that this manner of determining the charm coefficient 
function does not reflect the true physics, i.e. that a real charm-anticharm pair must be 
generated via the photon scattering, leading to the physical threshold in W'^, and cancellation 
between terms is required to reflect this correct threshold. A full critique of this approach can 



be found in |^6|, but we note here that our variable flavour number scheme is certainly very 
different to this alternative prescription (which is not yet implemented in their definition of 
NLO). 

The theoretical treatment of the bottom quark is essentially identical to that for the charm 
quark outlined above, and we use the FFNS below = ml and the VFNS above = ml. 



As discussed in [^, the procedure also generalizes to other processes in a simple manner. 

5.2 Implications for F| and the global analysis 

In Fig. ^ we show the result of the NLO calculation of F2 for x = 0.05. The very smooth 
transition from the description at low in terms of the FFNS to high in terms of the 
massless prescription demonstrates the success of our procedure. The result is qualitatively 
similar for all other x values. However, the obvious discrepancy between the continuation of 
the FFNS curve and the zero-mass curve for relatively low (i.e. ~ 5—20 GeV^), diminishes as 



^^We are grateful to Steve Riemersma and Jack Smith for providing the program to compute the 0(0^) 
contributions. 

^"We note that the same definition of the zeroth order coefficient function is adopted in although of 
course there are differences between this and [M, notably the mass dependent evolution in the former. 
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we go to lower x as we are then further from the physical threshold in W"^, and mass-dependent 
effects become less important. 

As was demonstrated in the global fit to structure function data achieved using our 
prescription for charm was superior to that using either the FFNS or the massless prescription. 
The former has too slow an evolution due to the lack of ln{Q'^/m'^) terms, and the latter gives a 
definite kink in F2{x, Q"^) at mass thresholds. Thus our global fit incorporates the best available 
treatment of heavy quark mass effects using NLO-in-ag QCD. 

One point to note is the influence of the charm prescription on the optimum value of as'(Ml). 
In the most recent global analysis this value came out to be 0.113, although an alternative set 
of partons was also given for ^^(Ml) = 0.120 p6[. The fact that the value of a5(M|) for 



the best fit has risen to 0.1175 is partially due to some of the new data in this fit, e.g. the 
reanalysed CCFR data, and the final NMC data. However, a fit to the new data using the old 
massless charm prescription results in a value of q;s(M|) = 0.116. Hence, the effect of the new 
treatment of heavy quarks is to increase the value of as{M^) by about 0.002. We also note 
that the prediction for Fl{x,Q'^) is very different to previous analyses, i.e. smaller; the charm 
contribution being very suppressed near threshold since the gluon coefficient function behaves 



like V . This leads to a better description of the NMC data on R [|1^, as seen in Fig. |18 



We can also look at the charm structure function directly. In Fig. ^we show the comparison 
of the charm structure function -^2^(0:, Q^) resulting from the MRST partons with all available 
data. The data at intermediate x values come from EMC measurements of inclusive muons, 
while the new data from HERA are obtained by measuring D and D* cross sections [^, ^. We 
show MRST curves together with those resulting from the MRST ((7 |) and MRST ((7 [) parton 
sets. In all cases the value of is taken to be 1.35 GeV. (A comparison for different values of 
rric can be found in |^B|.) The predictions using the central gluon agree well with the data. As 
we would expect the difference between the predictions is only significant for relatively high x. 
However, we see that for the two highest x bins the curves for MRST((yf |) tend to fall below 
the data. Thus, the EMC charm data are capable of acting as a weak constraint on the form of 
the gluon at high x, and it is clear that any parametrization which has an even smaller gluon 
at high X than MRST ((7 J.) would be inconsistent with these charm data. 

We should also justify our choice of tjIc = 1.35 GeV. This value is chosen somewhat as a 
compromise. As shown in |^ the charm data on their own prefer a value in the region of 
1.5 GeV. However, the quality of the global fits obtained are slightly sensitive to the value of mc 
and we find that rric ~ 1.2 GeV, or even lower, actually gives the best fit. Practically all this 
sensitivity to rric comes from the NMC F2 data with x < 0.1. As already discussed, these data 
grow with more quickly than the theory predicts. If rric is lowered then the charm evolution 
is slowed down less by mass effects, and the slopes of the theory curves increase. This is seen 
clearly in Figs. |^ and |15| which demonstrate the improvement in the description, although it 
is clear that a problem remains. As already mentioned there are other effects which could be 
responsible for the apparent discrepancy with the observed value of dF2/dlnQ^ in this region, 
and hence we do not take this as strong evidence for a low value of rric- Indeed, as seen in 
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|[46|| , the F2 data themselves completely rule out such a low value of nic. Hence we choose 
rric = 1.35 GeV as a compromise between the values required by the best global fit and the 
best description of charm data. 

Of course there are as yet no data on the bottom quark contribution to the structure func- 
tion, and because it contributes with a charge squared of 1/9, and only at relatively high Q"^, 
it forms only a very small fraction of the total structure function. Hence, the value of rrih has 
essentially no impact on the quality of the fit, and is always taken to be 4.3 GeV. Future data on 
the bottom structure function could act both as a verification of our treatment of heavy flavours 
and as a determination of the bottom quark mass. The b quark distribution at = 10^ GeV 
is shown in Fig. ^. 



6. Drell-Yan production 

The observation of Drell-Yan production in high energy pN collisions offers a valuable con- 
straint on the sea quark distributions since the leading-order subprocess is qq ^ 'j* ^ 
The small xp data of the E605 collaboration |jl5| are used in the global analysis and they 
constrain the sea in the interval 0.15 < x < 0.4. The description of the E605 data is shown 



in Fig. |28|. The curves are obtained using NLO QCD, with the factorization and renormal- 
ization scales set equal to the invariant mass M of the lepton pair, together with an overall 
phenomenological normalization parameter which allows for possible higher-order effects. 



The more recent measurements of the E772 collaboration span a larger kinematic range 
and, in principle, allow the sea quark distribution to be probed down to x ~ 0.025. In general 
the agreement between the data and the MRST prediction is quite satisfactory, see Fig. P^ . 
However, we point out a discrepancy at high xp and low ^/T. In this region the dominant 
contribution is 

da ~ u{xi) u{x2) + d{x2) (32) 

where Xi ~ xp and X2 — t/xi, and where the partons are sampled at scales = = st 
with ^/s ~ 40 GeV. If we assume that u{xi) is known, then there is a factor of two discrepancy 
between the E772 data and {u + d) evaluated at X2 ~ 0.03. We find that there is no way in 
which the global fit can remove this discrepancy since, at such small x values, {u + d) gives 
approximately a third of the total contribution to F2, and so is well determined by the deep 
inelastic scattering data. 



7. Flavour decomposition of the sea 

As seen in Figs. ^ and ^ the sea quark distributions (m, d, s, c and h) have an interesting 
non-trivial structure. Of course, as increases, all these distributions ultimately evolve to a 
common form, concentrated at small values of x, since they are driven hj g ^ qq transitions. 
However, at accessible values the different flavour sea quark distributions are quite distinct. 
Due to their heavy mass, the c and b can be treated perturbatively, as discussed in Section 5. 
The X and dependence of these heavy quark densities are therefore completely determined. 
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with their mass being the only free parameter. On the other hand for the hght quark distri- 
butions, u, d and s, we may use massless evolution, but here the distributions at the starting 
scale have a non-perturbative origin and are determined by experiment. For the light quark 
distributions we use the parametric forms given in (^, (^-(^. 

7.1 The strange quark distribution 

The observations of deep inelastic dimuon production indicate that s has the same x shape 
as {u + d) but with an overall suppression of the order of 50% at — 4 GeV^, presumably 
due to the mass of the strange quark. We reflect this behaviour by using an s parametrization 
which has the same shape as {u + d) at Qo, but with a single overall parameter (0.2 in (^) 
which is adjusted to fit to the CCFR dimuon data [0. Fig. |30| shows the MRST strange 



quark distribution compared to the measurements determined by the CCFR collaboration in a 



NLO analysis |16| of their dimuon production data. We see that our minimal parametrization 
gives excellent agreement with their experimental result, which supports the assumption that 
s/ {u + d) is essentially independent of x. As noted earlier, at small x, the differences between 
u^d and s decrease with increasing Q^, due to the dominance of the g ^ qq subprocesses. 
Finally we note that an independent measurement of s could, in principle, be made at HERA 
using the charged current subprocess e~ + s ^ + (c — > /i^) provided that the accelerator 
integrated luminosity was sufficiently high. 

7.2 Determination of the difference u — d 

The structure function measurements (of F^*^, F2"', and xF^^) determine (m + J), but 
not {u — d). The sum rules do give some information on the integral over u — d, which indicate 
that, on average, d is greater than u. 



For a direct determination of u — d we must look elsewhere. One proposal is to measure 
the asymmetry of Drell-Yan production in pp and pn collisions 

A — ^pp ~ '^P" \ — r 

^DY = 7 = (.>5'3j 

where r = (Tpn/cTpp and where a = d'^a/dMdxp with M and xp being the invariant mass and 
the Feynman x of the produced lepton pair. At leading order we have 

^ _ C^pn _ {^Uid2 + diU2 + ^Uld2 + dlU2 + 23132 + 8C1C2) 

app (4-U1-U2 + did2 + 4-uiM2 + did2 + 2siS2 + 8C1C2) 
where the 1,2 subscripts indicate that the partons are to be evaluated at 



XuX2 2 



^ i^±XF + ^xl + 4r) , (35) 



with T = M'^/3. We may rearrange the expression for 1 — r, and hence that for A-^y, to show 
that it is dependent on the combinations {ui — di) and {u2 — (^2)- 



20 



The first experiment of tliis type was performed by the NA51 collaboration [11]. They 
measured 

Rdp^^ = -il+r) (36) 

at xi = X2 = 0.18 and found Ady = —0.09 ± 0.02 ± 0.025, which corresponds to rf/w ~ 2. 
Very recently the E866 collaboration [|I^ have measured Rdp over a much wider range of M 
and xp, which enables a study of the x dependence of {u — d) over the range 0.04 < x < 0.3. 
The continuous curve in Fig. |3l| shows our fit to these dataQ The dotted curve shows the 
values which would have been obtained for the ratio if we were to set u equal to d, that is if 
we were to take = in (|^). The implications for d/u from our fit to the E866 data are 
shown in Fig. ^ Interestingly the structure of d/u shows that, at the maximum value of x 
that is measured^, the ratio has decreased to give d u. Moreover we see that the NA51 
measurement occurs at a value of x for which d/u is essentially at a maximum. Nevertheless 
the new data indicate a somewhat smaller value of d/u at this point, x = 0.18. For comparison 
we also show the prediction for d/u obtained from the MRS(R2) set of partons |^ — partons 



which were obtained from a global fit which included the NA51 measurement, but for which 
the E866 data were not available. 

Independent information on the u,d flavour asymmetry is currently being obtained by the 



HERMES experiment at HERA from the observation of semi-inclusive deep inelastic events. 
By observing final state vr^ mesons, they measure the ratio 

^ ajep ^ e-K-X) - ajen ^ e-R-X) 
^ ' ' a{ep^en+X) - a{en-^en+X) ^ ^ 

where z is the fractional energy of the fragmenting parton that is carried by the pion. The 
HERMES semi-inclusive data lie in the kinematic range 0.02 < a; < 0.3 and 0.3 < < 
10 GeV^. At leading order r{x, z) is a direct measure of {u — d)/{u — d) since 



l + r(x,z) {u-d) + {u-d) 2, + F 



1 — r(a;, z) (u — d) — (u — d) 5 \1 — F 



(38) 



where F{z) is the ratio of the disfavoured to favoured u ^ n fragmentation functions, F = 
Z)^ / D'^^ . Fig. ^ compares^ the MRST predictions with the preliminary HERMES measure- 
ments of {d — u)/{u — d) as a function of x. The good agreement between the MRST curve 
and this independent measure of m — c? is confirmation that d/u is now reliably known as a 
function of x and for x < 0.3. 

^^We calculate the ratio Rdp from the NLO QCD expression for cPa/dxidx2, with xi computed from the 
mean value of X2 for each bin, and then we join the resulting values of Rdp to form a smooth curve. 

^^Our parametric form should not be extrapolated to predict d/u at larger values of x, where both u and d 
are very smalL 

^^Note that we compare NLO MRST partons with a LO ratio extracted from the semi-inclusive data. However 
the effect of the NLO corrections is expected to largely cancel in the ratio. 
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Historically the first indication of the u ^ d flavour asymmetry of the sea came from the 
evaluation of the Gottfried sum 



/■I dx ,^ 
Jo X 



'GS 



FD 



(39) 



by NMC [51]. The final NMC measurements |52| of {F2^ — F2^) are compared with the values 
obtained from the MRST partons in Fig. 53. The area under the MRST curve yields the value 
0.266. This is slightly larger than the value Iqs = 0.235 ± 0.026 found at = 4 QeV^ 
The small discrepancy is induced in part by the requirement that the MRST 



'GS 



by NMC [H 



partons also fit the E866 data. 



8. W rapidity asymmetry 

The charge asymmetry at the Fermilab pp collider, 

daiW+)/dy-daiW-)/dy 
^^^^ " da{W+)/dy + da{W-)/dy' ^ 

is a sensitive probe of the difference between u and d quarks in the x ~ 0.1, Q ~ M]y region 
53| . Because the u quarks carry more momentum on average than the d quarks, the bosons 



tend to follow the direction of the incoming proton and the W~ bosons that of the antiproton, 
i.e. Aw > for y > 0. Thus a precise measurement of the W asymmetry serves as a valuable 
independent check on the u- and d-quark distributions. In practice it is the lepton asymmetry, 

^^^^^ = ain + aie-y 

which is measured, where cr{i^) = da/dy^ are the differential pp —>■ W^X i^uX cross 
sections for producing leptons of rapidity yf . There is a direct correlation between the 
lepton asymmetry and the slope of the d/u ratio. To see this we first note that the dominant 
contribution to production comes from the ud {du) annihilation process. Thus 

_ u{xi)d{x2) - d{xi)u{x2) .^2) 
u{xi)d{x2) + d{xi)u{x2) ' 

where the scale Q = Mw is implicit for the parton distributions, and 

2^1,2 = Xoexp{±y), xo = (43) 



If we introduce the ratio Rdu{x) = d{x)/u{x), then, for small y, 

Aw{y) ^ ~xoy (44) 

where the prime denotes differentiation. In reality, the situation is of course more complicated 
— it is the lepton asymmetry which is measured, and there are subleading and higher-order 
corrections to (0). Nevertheless, the correlation implied by p^) is evident in the full prediction. 
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The CDF collaboration [|T2] have recently extended the range and improved the precision 
of their measurements of the asymmetry Ag. The new data extend to higher values of lepton 
rapidity yi and the measured values are below the extrapolation of the previous global fits. 
The fitQ to the lepton asymmetry data is shown in Fig. For comparison we also show the 
result from a previous set of partons, MRS(R2), which were fitted to earlier CDF asymmetry 



measurements. Recently the effects of soft gluon resummation on Ai have been calculated ||55 
They increase the asymmetry slightly at the highest yi values of the data. For example at the 
highest value, ye = 2.2, the increase in A^ is about 0.02 so that the MRST prediction would, 
as it happens, be raised to coincide exactly with the data point. 

In order to accommodate the new Ai measurements the d distribution increases with respect 
to the M at a; ~ 0.3 (in comparison with our previous global analysis). The change affects 
much more than Ff'. The consequence is that the ratio F2 /F2 is increased, and the description 



of the new NMC data is improved relative to MRS(R2), see Fig. ^ 



9. Implications for hadron collider cross sections 

According to the QCD factorization theorem, 'hard scattering' hadron collider cross sec- 
tions can be expressed in terms of parton distribution functions convoluted with perturbatively 
calculable subprocess cross sections, 

dax = J2 jdxidx2fi{xi,^^)fj{x2,^i^)ddij^x, (45) 

with, for example, X = W^, Z,QQ, jets or Higgs. Processes with well-measured final states, 
and for which the next-to-leading order corrections to a are known, can therefore provide impor- 
tant cross-checks on the parton distributions. A prime example is the W rapidity asymmetry 
at the Tevatron which, as we have just seen, is used in our global analysis to constrain the u 
and d distributions. Leading-order kinematics imply xi,2 = Mx exp{±yx) / \/s. In general, at 
colliders like the Tevatron and the LHC, quark and gluon distributions are probed at x values 
where they are 'measured' by deep inelastic scattering and prompt photon experiments, but at 
significantly higher scales fi"^ ~ M]^. 

A precise knowledge of parton distributions is absolutely vital for reliable predictions for 
signal and background cross sections at the LHC [^]. Uncertainties can arise both from the 
starting distributions and from DGLAP evolution. A detailed assessment of these uncertainties 



illustrated with reference to various standard cross sections will be presented elsewhere [57]. In 



this section we focus on several standard cross sections measured at the Tevatron pp collider. 
9.1 W, Z production 

While the W (lepton) rapidity asymmetry probes the relative size of the u and d distribu- 
tions, the total cross sections for W and Z boson production in pp collisions at ^/s = 1.8 TeV 

^^The curves in Fig. ^ are calculated using the next-to-leading-order program DYRAD of ||5^ . We thank 
Nigel Glover for helping with these calculations of the W asymmetry. 
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provide an important check of the overall magnitude of the quark distributions in a region of 
X ~ 0.05 — 0.1 where they are constrained at lower /i^ by deep-inelastic (in particular NMC) 
data. Since the perturbative QCD subprocess cross section is known to next-to-next-to-leading 
order |^ and the electroweak parameters are precisely determined, there is very little theoret- 
ical uncertainty in the predictions once the parton distributions are specified. 

We begin by displaying in Fig. the W and Z total cross sections times the leptonic 
branching ratios^ as measured by UAl and UA2 |Q at -y/s = 630 GeV and by CDF 



and DO ||62| at y/s = 1.8 TeV, together with the predictions obtained using our default MRST 
set.Q Only the most recent (Run lA) Tevatron published measurements are included. The 
factorization and renormalization scales are set equal to and Mz respectively. The overall 
agreement between theory and experiment is excellent. 

To study this in more detail, we focus on the more precise Tevatron measurements and 
show, in Fig. |3^, the predictions of the five MRST sets (the default MRST set, together with 
(7 1,(7 I, as TT)Ct5 ii)- There is an overall spread of approximately ±2% about the default 
prediction,^ significantly smaller than the current experimental errors. The variations in the 
predictions are easily understood. At these (small) x values, the quark distributions increase 
with increasing /i^. The larger the 0:5 the faster the increase, and so <7v{(ys ii) < crv'(MRST) < 
evicts TT)- The different gluon distributions also give rise to differences in the ay predictions. 
In this X ~ 0.05 region, the ordering of the gluon distributions is (7 j < (^(MRST) < g I, 
see Fig. The larger the gluon the more rapid the DGLAP evolution, and so avio T) < 
(Tv'(MRST) < cry ((7 i). The gluon variation is slightly smaller than the as variation. 



We may conclude from Fig. |3^ that for a given set of quark distributions fixed by DIS 
data at a relatively low ~ 10 GeV^ scale, the net uncertainty in the aw,z predictions at 
high ~ 10^ GeV^ coming from the gluon, as and unknown higher-order corrections is of 
order ±2%. The true theoretical uncertainty has a significant additional component from the 
uncertainty in the absolute normalizations of the {u and d) distributions as determined by 
the normalization uncertainties in the structure function data themselves. In the relevant x 
range for the Tevatron W and Z cross sections, the main constraints on the quarks come from 
the precise NMC and F2 data. These have an overall systematic (normalization) error of 
approximately ±2.5%, which would give a corresponding error of ±5% on the weak boson cross 
sections, comparable to the current experimental errors from CDF and DO. However this is 
almost certainly a large overestimate since the requirement of mutual consistency between the 
various deep inelastic data sets, together with the sum rule constraints, does not permit the 
possibility of changing the normalization of an individual data set over its full allowed range. 



This point will be addressed in ||57[| . 



25The Standard Model values B{W Iv) = 0.1084 and B{Z 1+1-) = 0.03364 are used. The electroweak 
boson masses are taken to be Mw = 80.43 GeV and Mz = 91.1887 GeV. 

Strictly speaking there is a slight inconsistency from combining a subprocess cross section calculated to 
NNLO with NLO-evolved partons. However the NNLO contribution is numerically very small. 

^"^A measure of the scale dependence of these predictions is obtained by using instead n = My 12 and 2My. 
The effect is shown as error bars on the default prediction and is evidently very small. 
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A final point concerns the impact on the experimental cross sections of the luminosity 
measurement and uncertainty. The latter (±3.6% for CDF and ±5.4% for DO |^T], ^) dominate 



the most precise {W eu) cross section errors. In addition, the value assumed for the total 
pp cross section is slightly different for the two experiments, and this may account in part for 



the systematically smaller DO cross sections displayed in Fig. |38 



9.2 Top production 

At the Tevatron, ti production occurs dominantly (~90%) via qq annihilation. The cross 
section therefore samples valence u and d quarks at x > 2mtl \fs ~ 0.2 — 0.3 and ^ ~ 



10^"^ GeV^. The leading order subprocess cross section is proportional to a|, but the en- 
hancement obtained by increasing 05 is partially compensated by the decrease in the parton 
distributions at large x caused by more rapid DGLAP evolution. 

The production cross section is known exactly at NLO [Q. It is traditional to estimate the 
residual theoretical scale dependent uncertainty by varying the factorization and renormaliza- 
tion scales in the range mijl < /i < Ivnt which, at the Tevatron, gives an approximate ±10% 
variation about the [i = vrit prediction for a fixed set of partons and fixed as- 

Another important effect is the higher-order contributions from soft gluon emission, which 
are expected to be large when the tt system is produced near threshold. Techniques have 
recently been developed for resumming the dominant leading (LL) ln(l — 4mVs) logarithms 



pl |. In Ref. [^] the NLL logarithms have been computed and used to obtain an 'improved' 
resummed cross section. The effects of resummation are particularly large for large scale choices, 
whereas for /i = mt/2 the NLO cross section approximates the resummed cross section to better 
than 1% at Tevatron energies ||65||p^. In our calculations we will therefore compute the top cross 



sections using NLO QCD, with mf = 175 GeV and /i = mt/2. 



Figure |3^ shows predictions for the total tt cross section at the Tevatron using the five 
canonical MRST parton sets. Data points from CDF and DO |^ are also shown. The cross 
sections are slightly larger for the sets with the larger (large-x) gluon and larger as (although 
for the latter, note that the variation is smaller than the naive estimate of ±10% from the 
change in the overall q;| value would suggest). The uncertainties due to the valence quarks, the 
gluons, and as are all at the ±5% level. The resulting overall parton distribution uncertainty 
in the (T(tt) prediction is therefore at the ±10% level, comparable in magnitude to the scale 
dependence (combined with other higher-order effects ||6^) uncertainty. The agreement with 
the current CDF and DO measurements is entirely satisfactory. 

9.3 Large jet production 

The single jet inclusive distribution at the Tevatron is a particularly interesting observ- 
able. Even though the measured cross section falls by more than six orders of magnitude when 
Et increases from 50 to 400 GeV, the NLO QCD predictions reproduce the data to well within 



^^For the choice [i — nit the effects of beyond-NLO resummation increase the NLO cross section by approxi- 
mately 5% ia. 
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the systematic error band. Nevertheless the detailed shape of the spectrum, taking only the 
statistical errors into account, shows interesting features. The spectrum has been measured by 
both the CDF and DO collaborations |5S[. The experimental spectra are compared with 



the NLO predictions]^ obtained from five of our canonical sets of partons in Figs. |D-^. In 
each case we use the prediction of the MRST parton set as the base line for the comparison. 
The predictions contain an overall normalization factor which is adjusted to give the optimum 
description of the data. The value of the factor is shown on the plots for each set of partons. 



Figs. ^ and ^ compare the CDF and DO spectra with the predictions obtained from the 
three parton sets based on different gluons. It is interesting to see that the set, MRST(5f |), 
with the larger gluon at large x (that is the set in which the WA70 prompt photon data are 
fitted with partonic (kx) = 0) gives the best description of the shape of the observed spectrum 
and an overall normalization nearest to unity. The effect is particularly pronounced for the DO 
data. Figs. ^ and ^ compare the CDF and DO spectra with the predictions of parton sets 
corresponding to three different values of as^M^). In this case the parton set, MRS(a5' J.|), 
with the smallest as gives the best description. 

Note, however, that the comparisons between experiment and theory shown in Figs. HO-EH 



should not be taken too literally. Aside from the large experimental systematic errors which 
have not been included, there are residual theoretical uncertainties, for example from scale 
dependence and the precise modelling of the experimental jet algorithm, which are important 
at the (9(±5%) level (see for example the recent study in Ref. ||6^). It is therefore premature 
to draw any definite conclusions about preferred gluons distributions and/or as values. 

A final point concerns the impact of 'intrinsic' transverse momentum smearing on the jet 
Et distribution. Using the same procedure as implemented in the prompt photon studies of 
Section 4, we have investigated the effect of different choices of {kx) on the shape of the distri- 
bution. Qualitatively, the effect is to steepen the distribution slightly at the low Et end of the 



spectrum, as for the Pt{i) distributions shown in Figs. 23-25. For example, for {kx) = 4 GeV 
we find that the jet cross section is increased by +3% for Et = 75 GeV, and by less than +1% 
for Et > 150 GeV. The description of the CDF and DO data by the canonical MRST parton 
set is not therefore improved. 



10. Conclusions 

In order to examine the partonic structure of nucleons we have performed a global analysis 
of data on deep inelastic scattering and related hard scattering processes using NLO-in-a^ 
QCD. This present treatment represents a significant improvement over our previous analyses 
for a number of reasons. First, there is the availability of a number of updated, or completely 
new sets of data. These have all been discussed in the introduction. However, we note that 
those experiments which have had a major impact on the changes to the parton distributions 
include those that produce a relatively small number of data points: i.e. the new data from 

We are grateful to Nigel Glover for help in performing these calculations. 
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the E866 and HERMES experiments, which constrain the value of {u — d); the extended data 
on the W rapidity asymmetry, which constrain d/u; and the new prompt photon data from the 
E706 experiment, which probe the gluon at high values of x. 

There are also new features in the way in which we perform the analysis itself. For the first 
time we incorporate intrinsic parton when examining those data which are sensitive to it, 
which in practice means the prompt photon data. We also use a new, much more sophisticated 
treatment of the charm (and bottom) contribution to structure functions. This new procedure 
naturally includes both a smooth behaviour in the threshold region and the summation of large 
logarithms at high Q^, and gives a parameter free (up to the charm mass) prediction for the 
charm contribution. For the first time we have also examined the sensitivity of the partons sets 
obtained to cuts in both and x imposed on the data, letting the former vary between 2 GeV^ 
and 10 GeV^ and the latter up to x^i^ = 0.01. We have also investigated the uncertainties 
which exist for our determination of the strong coupling constant 05 in a more systematic 
manner than in previous analyses, obtaining a quantitative estimate of the allowed variation 
about our central value of ^^(Ml) = 0.1175. Finally, for the first time we have also made a 
thorough investigation of the uncertainty in the gluon distribution obtained from our analysis, 
producing two sets of partons which represent two extremes as well as the central, preferred 
set. 



Both the new data sets and the new theoretical procedures have led to significant differences 
between the MRST partons and those produced by previous analyses. In order to exhibit these 
differences we can look at Fig. ^ which shows the comparison between the MRST partons 
and the preferred set from our previous global analysis, i.e. MRS(R2). First, looking at the 
comparison at = 20 GeV^ (upper plot), most partons show significant differences. The new 
gluon is much smaller for x ^ 0.2, since the inclusion of intrinsic means a smaller gluon is 
required to fit the WA70 prompt photon data. The form of the gluon at small x is constrained by 
the HERA data, and therefore the momentum sum rule allows a larger gluon at x ~ 0.05. The 
form of the charm quark distribution is clearly very different to our previous MRS(R) analysis 
26| , where the evolution of the charm quark in the MS scheme took place from Q'^ = 1 GeV^ 



rather than the correct value of Q = m^, and a phenomenological damping factor was used. It 
is this damping factor which led to the previous small charm distribution at small x. However, 
the charm distribution is driven by the gluon distribution, and one can see that the difference 
in shape between the new and old charm distribution mirrors that of the gluon. It is this which 
causes the MRS(R2) charm quark to be larger at large x, overcoming the effect of the damping 
factor which is a function of only. The light quark distributions are all a little larger for 
the MRS(R2) partons at small x, compensating for the smaller charm distribution. At high 
X the u quark is essentially unchanged, the d quark is larger in the MRST partons in order 
to accommodate the extended W asymmetry data, and the strange quark (which effectively 
represents the sea quark distribution) is a little smaller, presumably due to the requirements 
of fitting the Drell-Yan data with a slightly larger valence quark distribution. Also examining 
the lower plot of Fig. ^ we see that the systematic differences between the light partons are of 
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the same form as at the lower Q"^ value, but have been washed out somewhat (in particular the 
small X form of the partons is almost identical in each set). However, the charm distribution 
is still very different. Again, the shape difference mirrors that of the gluon, but the MRS(R2) 
charm distribution is larger. This is because the damping factor is now unity at such a high 
value of Q^, but the fact that the evolution of the MRS(R2) distribution began at lower Q'^ 
results in a constant difference between the two distributions, i.e. the charm generated in the 
evolution from Q"^ = 1 GeV^ to Q'^ = m^. Although this difference disappears in the ratio 
asymptotically, it is still significant at = 10* GeV^. 

We can also compare the new MRST partons with the preferred set of the most recent 
CTEQ analysis, i.e. CTEQ4M. This comparison is displayed in Fig. We first compare 
the partons at = 10 GeV^. In the high x region there is very little difference between 
the CTEQ4M partons and the MRS(R2) partons. The CTEQ4M strange (and therefore sea) 
distribution becomes much larger at x > 0.4, but this is beyond the range of the E605 Drell-Yan 
data. Hence, all differences between quarks in this region are presumably due to the new data 
included in the MRST analysis. Even the gluon is similar to MRS(R2) in this range, being 
constrained mainly by the high-i?T jet data, which, as we have already noted, require a high x 
gluon similar to MRST((7 |). However, the CTEQ4M gluon becomes significantly larger than 
the MRST gluon at small x, leading to a stronger growth of F2{x,Q'^) with in this region. 
This is partly due to the cut of 4 GeV^ imposed by CTEQ, and a similar trend is seen for 
our gluon when similar cuts are imposed, as demonstrated in Fig. Since the evolution of 
the charm parton distribution is the same for CTEQ4M as for MRST (the treatment of the 
coefficient function being very different), differences in the charm are largely due to differences 
in the gluon, and indeed the shape of the charm ratios mirrors that of the gluon ratios. The 
light quark distributions at small x are a little smaller than those for the MRST partons, but 
as with the gluon differences this can be attributed mainly to the different cuts, the same 
qualitative effect again being seen in Fig. [2^. However, the slightly different ratio between the 
u and d quarks and the strange quarks, even at x = 10~^, is due to the CTEQ requirement that 
2s/ {u + d) = 0.5 at = 2.56 GeV^ rather than 1 GeV^, which results in the strange quark 
being a little smaller. In the intermediate x region the CTEQ4M u and d quarks are both 
slightly larger than those for MRST, compensating for the smaller strange and charm quarks. 
As with the previous comparison, differences between the parton sets tend to be washed out at 
very high Q^. 

In summary the main differences between the CTEQ4M and MRST partons arise from (i) 
our inclusion of new data, (ii) the different treatment of the gluon at large x and (iii) the 
choice of cuts in Q^. However, there are some small additional systematic differences between 
the two analyses. The normalizations imposed on the data sets are a percent or two higher 
for CTEQ than MRST. If we compare CTEQ4M partons with those obtained in our analysis 
using a comparable cut (Q^ > 5 GeV^) then the partons are more similar, but significant 
differences remain. These may be due to the systematic difference between the output of the 
evolution programs, as noted in |69]. With the quality of data now available such differences 
are becoming significant. 
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Let us briefly justify our confidence in our new canonical set of partons (MRST). As we 
have seen they give an excellent overall description of the data. There are certain conflicts: 
as in our previous analyses we omit the (new) CCFR neutrino F2 and xF^ data for x < 0.1 
from our fit as explained in Section 3; and we have also noted a discrepancy with the E772 
Drell-Yan data for large xp and small r. There is also a slight systematic discrepancy with 
the measured value of (x, Q^) / dln.Q'^ for the NMC data below a; = 0.1 and for some of the 
HERA data. The quality of the fit to these data is clearly very satisfactory, but this observed 
effect may be a sign that NL0-in-a5 DGLAP evolution is not sufficient at small x, and that 
theoretical corrections are required.^ It will be interesting to see whether further HERA data 
confirms the observed trend, and measurements at high y may be particularly important since 
they probe both F2(x,Q^) and Fl{x,Q'^). We see from Fig. |21] that if we vary as away from 
our central value of 0.1175 in either direction, then the quality of the fit to one or more of the 
data sets becomes rather poor very quickly. Hence, variation of allows very little variation 
in our partons. We believe that our default set of partons, denoted simply by MRST, gives the 
most reliable treatment of the gluon at high x. The E706 prompt photon data have made it 
clear that intrinsic parton kj- is required, as one might expect from theoretical arguments and 
from the evidence already observed in Drell-Yan data. It is also very plausible that the amount 
of intrinsic kp increases with centre of mass energy, and this hypothesis has recently received 
support from an approximate soft-gluon resummation calculation In the absence of a fully 
rigorous theoretical treatment of the broadening of the prompt photon px distribution we have 
adopted an experimentally motivated Gaussian form for the parton intrinsic kx smearing. The 
MRST(5f I) gluon obtained from assuming zero intrinsic fcy for the fit to WA70 data is clearly 
at one extreme, while the MRST(5f |) gluon obtained from allowing the partons to have the 
maximum {kx) in the fit to WA70 data is at the other extreme by construction. The default 
MRST gluon has (fcr) intermediate to these extremes. While we cannot be absolutely sure of 
the most appropriate values of {kj) to take for the two prompt photon experiments they cannot 
be much different to those employed in MRST. Since MRST((yf |) and MRST((yf [) represent 
such extreme variation in parton (/ct) we feel that the variation in the gluon they encompass is 
also sufficient to account for the smaller effect of scale dependence in prompt photon production. 
We note that there is no dependence on intrinsic parton kx in any of the other data sets we 
analyse except for the high Ft jet data from CDF and DO, where the effect is very small, and 
not helpful. 

We now have the most complete determination of the proton parton distributions yet 
achieved, both in terms of the theoretical treatment and in terms of the data analysed. How- 
ever, in both there are possible improvements to be made. In terms of the theoretical procedure 
used in the analysis there are two major points. At present we use a leading twist NLO-in-a^ 
calculation. However, there are both leading ln(l/x) and higher twist corrections which should 
have a significant effect in certain kinematic domains (we impose a cut in W"^ in order to re- 
move the sensitivity to higher twist effects at high x and low Q^), and a hint of this is perhaps 

•^"The discrepancy with the E772 data occurs when one of the partons has x ~ 0.03, but it seems unhkely 
that such a large discrepancy can be cured by leading ln(l/a;) effects at such a high value of x. 
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seen in the results of our analysis when different cuts are applied. In order to be truly confi- 
dent of our results over the whole kinematic domain we need either proof that these possible 
corrections are very small, or a reliable method of including them in our analysis. Also, our 
present treatment of intrinsic parton is phenomenologically inspired. While we believe it 
is a good representation of the true effects we await an improved theoretical understanding of 
the origin of both perturbative and nonperturbative intrinsic kx, hopefully pinning down the 
actual amount in a quantitative manner. 

There are a number of improvements in data which would also be very important in in- 
creasing the accuracy of determinations of parton distributions. Obviously, any increase in the 
precision of existing data and/or a widening of the kinematic range would be useful. It would 
be advantageous to have further sources of neutrino structure function data as a consistency 
check, since at present the analysis relies on a single high precision experiment. Of particular 
use would be an increase in both the amount and precision of data on the charm structure func- 
tion, as this would help differentiate between different prescriptions for heavy flavour effects, 
be a direct probe of the gluon, and possibly be important in determining small x corrections. 
We anticipate data of this type from HERA in the near future. Likewise, data on the bottom 
structure function would be useful for the same reasons, but since the mass is higher, it would 
be a cleaner probe of perturbative QCD. Finally, it would be very interesting to have a direct 
measurement of Fi{x,Q'^) at small x, as this is an excellent test of different approaches to 
incorporating small x corrections. 

In summary we have used all the deep inelastic and hard scattering data in order to obtain 
what we believe is the most precise determination of parton distributions of the proton to date. 
On the theoretical side this illuminates the short-distance structure of the proton, while at 
the same time it provides an essential ingredient for extracting new physics from current and 
forthcoming hadron colliders. 
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TABLE I. The numerical values of the parameters of the starting distributions of three parton 
sets which differ in the value of the initial state partonic (/ct) used to describe the prompt photon 
data. The three columns correspond to (/ct) = 0.64, 0.4 and GeV for WA70 respectively. 
Note that Ag is fixed by the momentum sum rule and that A^^A^ are fixed by flavour sum 
rules, and are not therefore free parameters. 
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TABLE II. The fractions of the total momentum of the proton carried by the various partons 
in the MRST set. 
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III. Processes studied in the global analysis (* 
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TABLE IV. The values for the DIS data included in the three global fits which resulted in 

the parameter values listed in Table I. 



Data set 


No. of 

data pts 
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MRST(^ t) 
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164 
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66 
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67 


NMC n/p 


163 
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174 



The relative normalizations of the data sets are taken to be unity, except that the 
BCDMS /xp data is normahzed down by 2% and the SLAC ep data up by 2.5%. 
The CCFR data are corrected for heavy (iron) target effects using the information 
obtained from the muon-nucleus measurements. Since the a; < 0.1 corrected CCFR 
data are in disagreement with the NMC data, they are omitted from the fit. Only 
statistical errors together with a 1.5% uncertainty (to represent uncertainty in the 
heavy target correction) are included in the for the CCFR data since no overall 
systematic errors are given. 

All deuterium data are corrected for shadowing effects using the method of Badelek 
and Kwiecinski p9[ . 



Notes: (i) 
(ii) 

(iii) 
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Figure 1: The gluon distributions at = 2, 5, 20 and 100 GeV^ corresponding to the MRST, 
MRST(gi I) and MRST((/ |) sets of partons with, respectively, the central, larger and smaller 
gluon at large x. We take MRST as the default set of partons throughout the paper. 
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Figure 2: The ratios of the g ] and g [ gluons to the 'central' gluon (MRST) at = 10 and 
10^ GeV^. For comparison the ratio of the CTEQ4M |^T| gluon to our central gluon is also 
shown. 
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Figure 3: MRST partons at ^ 20 GeV^. 
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42 




Figure 5: Comparison of the MRST partons with those of the previous MRS(R2) set at Q'^ — 
20 GeV^ (upper plot) and ^ GeV^ (lower plot). 
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Figure 6: The description of the Ff" data at small x by the MRST set of partons. The 
comparison is made at twelve values of x chosen to be the most appropriate for the new 
HERA data. For display purposes we add 0.5(12 — i) to -Pf each time the value of x 
is decreased, where i = 1, 12. The experimental data are assigned to the x value which 
is closest to the experimental x bin. Thus the ZEUS data are shown in groupings 
with X values (3.5,6.3,6.5 x IQ-^), (1.02,1.20 x 10"^), (1.98,2.53 x 10"^), (4.0,4.5 x 10"^), 
(6.32,8.00 X 10-^), (1.02,1.20 x lO'^), (1.612 x 10-^), (2.53,2.60 x lO^^)^ (4 00 x lO^^), 
(6.325 X 10^3), (1.02 X lO'^), (1.612 x lO^^), and the HI data ^ in groupings with x val- 
ues (3.2,5.0 X 10"^), (0.80,1.30 x 10"^), (2.0,2.5 x 10"^), (3.2,5.0 x 10"^), (6.3,8.0 x 10'^), 
(1.3x10^3)^ (1.585x10-3), (2.0,2.5,3.2x10-3), (3.98,4.0,5.0x10^3), (6.3 xlO^^), (8.0x10^3), 

(1.3 X 10-^). The E665 data p4[, which are shown on the curves with the five largest x val- 
ues, are measured at x = (2.46 x 10-3), (3.698,5.2 x 10-3), (6.934 x 10-3), (8.933 x 10-3), 
(1.225,1.73 x 10-2). 
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Figure 7: Comparison of the F2 predictions of the MRST partons with the prehminary 1995 
and 1996 nominal vertex data of HI |]2^. For display purposes we add 0.2(19 — i) to F2 each 
time the value of x is decreased, where i = 1, 19. These data are not used in the global analysis. 



45 




10 10^ 10 10^ 



(GeV^) 

Figure 8: Description of the NMC F2 data by the MRST partons. The effect of lowering 
the charm quark mass from 1.35 to 1.20 GeV is shown by the dotted curve. For display purposes 
we have multiplied by the numbers shown in brackets. 
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Figure 9: Comparison of from the MRST set with the NMC deuteron data |T^. The effect 
of lowering the charm quark mass from 1.35 to 1.20 GeV is shown. For display purposes we 
have multiplied by the numbers shown in brackets. 
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Figure 10: The description of the large x BCDMS [23] and SLAC |25| measurements of by 
the MRST partons. 
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Figure 11: The description of the CCFR measurements of by the MRST partons. Only 
the data with x > 0.1 are included in the global fit. An x-dependent heavy target correction 
has been applied. The errors shown correspond to the quoted statistical errors together with 
a 1.5% uncertainty to represent the uncertainty of the heavy target correction. For display 
purposes we have added to the numbers shown in brackets. 
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Figure 12: The same as Fig. 11 but for the structure function xF^^ . 
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Figure 13: The effective exponents Xs, Xg of the small x behaviour of the sea-quark and gluon 
distributions of the default MRST partons versus Q^, defined such that xS ~ x~^^ and xg ~ 

X~'^^. 
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Figure 14: The description of the slopes dF2/dlnQ^ versus x for the HERA data. The experi- 
mental points were computed from linear fits to the data in Fig. The solid curve corresponds 
to the MRST set while the dashed line indicates the effect of removing low data (with 
< 10 GeV^) from the fit. The values of appropriate to the data are given (in GeV^) 
approximately by the formulae shown and the corresponding scale is indicated on the upper 
edge of each plot. 
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Figure 15: The continuous curve is the description of the slopes dF2/dlnQ^ of the NMC proton 
data. The experimental points were computed from linear fits to the data in Fig. |[ We also 
show the effect of modifying the charm quark mass and of excluding low x data (with x < 0.01) 
from the global fit. 
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Figure 16: The description of the quantity a = F2 — + (1 — compared with the 

prehminary HI 1995/96 ||2^ data. Also shown (dashed curve) is the contribution to a of F2 
only. 
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Figure 17: Prediction for Fl at HERA versus x for y = 0.6, 0.7, 0.8. Note that Fl has an 
imphcit y dependence via the relation Q"^ — xys. 
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Figure 18: Comparison of the predictions for R = Fl/2xFi with the NMC data. The values 
of appropriate to the data are given (in GeV^) approximately by the formula = 262x^'^^, 
and the corresponding scale is indicated on the upper edge of the plot. 



56 



X 



4 - 



Q^=2 GeV^ 

Fits using 
qS 2 GeV^ all x 
> 10 GeV^ all x 
qS 2 GeV^ X > 0.01 




10 



7.5 



5 - 



2.5 





20 




30 




Figure 19: The gluon distributions which result from making different cuts in Q'^ and x to the 
data included in the fits. The cuts which we use are specified in the first plot. 
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Figure 20: The parton distributions at = 10 GeV^, compared with the default MRST 
partons, obtained by making Qf — 5 GeV^ and Ql — 10 GeV^ cuts in Q'^ to the data included 
in the fits. 
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Figure 21: The contributions to the total global fit from the various data sets as a function 
of as- The parton set corresponding to the optimum value as — 0.1175 is denoted simply 
MRST and is the default set of partons used throughout the paper. The four other sets, which 
correspond to the adjacent values of as indicated by arrows, are used for comparison purposes. 
The values for the CCFR data are obtained from the statistical error and an additional 1.5% 
'systematic' error added in quadrature. 
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Figure 22: The extrapolation of the fits at x = 0.45 to high using the MRST, MRST(q;5 TT) 
and MRST (q;5 J, J,) set of partons. 
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Figure 23: Comparison of the WA70 [^] data with the MRST parton set with and without 
smearing in transverse momentum. 
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Figure 24: Comparison of the E706 data at 530 GeV with the MRST parton set. The 
results for two choices of the factorization scale are shown, together with the improvement 
of the description obtained by including parton transverse momentum in each case. The two 
'improved' curves lie on top of each other. 



62 




Figure 25: Comparison of the E706 data at 800 GeV with the MRST parton set. The scale 
is chosen to be pt/2 and the effect of including parton transverse momentum is shown. 
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Figure 26: The continuous curve is the MRST charm structure function ai x — 0.05 obtained 
using the NLO prescription and NLO evolution (VFNS). For comparison, the prescriptions of 
the fixed fiavour number scheme (FFNS) and the zero mass scheme (ZM-VFNS) are also shown. 
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Figure 27: The prediction for the charm structure function using the three MRST parton sets 
characterized by the different large x gluon behaviour. The HERA data at small x are from 
HI § and ZEUS §, and the large x data are from EMC §. 
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Figure 28: Hadroproduction of dileptons computed from the MRST parton set compared with 
the E605 data |T^. The theory curves inchide an additional K' factor of 0.9. No correction for 
the heavy target has been made. The scale on the left-hand axis is appropriate for the theory 
and data dX xp = —0.125. For display purposes the normalization is then decreased by a factor 
of ten for each step upwards va. xp- 
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Figure 29: Hadroproduction of dileptons computed from the MRST parton set compared with 
the E772 data [jl4[. The theory curves include an additional K' factor of 0.96. The scale on 
the left-hand axis is appropriate for the theory and data at xp = 0.025. For display purposes 
the normalization is then decreased by a factor of ten for each step upwards in x^. 
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Figure 30: Comparison of the strange quark distribution from the MRST set of partons com- 
pared with the determination of the strange sea obtained by the CCFR collaboration |T^ in a 
NLO analysis of their data on the neutrino production of dimuons. 
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Figure 31: The continuous curve is the MRST description of the E866 data for the ratio 
of the cross sections for hadroproduction of dileptons for proton and deuterium targets versus 
X2i the fractional momentum of the parton in the target. The other curves are for comparison 
only. 
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Figure 32: The ratio of the parton distributions d/u at = 30 GeV^ for the MRST and 
MRS(R2) parton sets compared to the estimate from the NA51 ||1T| measurement of the Drell- 
Yan asymmetry. 
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Figure 33: Predictions for the ratio {d — u)/{u — d) from the parton sets MRST and MRS(R2) 
compared with the prehminary estimates obtained by the HERMES collaboration from 
their measurements of semi-inclusive pion production in DIS. Note that there is an additional 
overall systematic error on the data of ±0.07. The values of the data vary with x approxi- 
mately according to the formula (GeV^) = ISx*^'^"^. 
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Figure 35: The description of the lepton asymmetry for production in pp coUisions at ^/s 
= 1.8 TeV. The data from CDF [0 are compared with the new MRST parton set and with 
the previous set MRS(R2). 
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Figure 36: The large x data on the ratio F^/Ff extracted from the measurements of F2 / F2 
by NMC |3[ compared with the MRST and MRS(R2) descriptions. 
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W, Z production cross sections 




Figure 37: Total W, Z production cross sections times leptonic branching ratios as a function 
of the (pp) colhder energy -^i, calculated using the default MRST partons. Experimental 
measurements from UAl UA2 |6^, CDF and DO ||6^ are also shown. 
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Figure 38: Predictions for the W and Z production cross sections times leptonic branching 
ratios in pp collisions at 1.8 TeV using the five MRST parton sets. The error bars on the 
default MRST prediction correspond to a scale variation of /i = My jl IMy, V = W,Z. 
Experimental measurements from CDF |^ and DO are shown. 



76 



10 



MRST partons 
NLO QCD, |x = m/2 

m. = 175 GeV 



CDF 



8 - 



DO 



7 - 



6 - 



5 - 



:10% 



4 - 



Figure 39: Predictions for the tt total cross section in pp collisions at 1.8 TeV using the five 
MRST parton sets. The error bars on the default MRST prediction correspond to a variation 
in the top mass of ±5 GeV. Experimental measurements from CDF ||6^ and DO are shown. 
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Figure 40: The next-to-leading order QCD description of the CDF single jet inclusive Et 
distribution by the MRST set of partons. The overall normalization of the QCD prediction is 
fitted to the data. The comparisons with the MRST ((7 f ) and MRST ((7 J.) parton sets are also 
shown. 
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Figure 41: The next-to-leading order QCD description of the DO single jet inclusive Et 
distribution by the MRST set of partons. The overall normalization of the QCD prediction is 
fitted to the data. The comparisons with the MRST ((7 f ) and MRST ((7 J.) parton sets are also 
shown. 
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Figure 42: As for Fig. |0 but including the comparisons with the MRST(a;5 It) and 
MRST(as ID parton sets. 
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Figure 43: As for Fig. il but including the comparisons with the MRST(a;5 It) and 
MRST(as ID parton sets. 
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Figure 44: Ratio of the partons of the CTEQ4M set to those of the MRST set at 
and 10^ GeV^. 
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